Smoothed universal correlations in the two-dimensional Anderson model 
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We report on calculations of smoothed spectral correlations in the two-dimensional Anderson model 
for weak disorder. As pointed out in (M. Wilkinson, J. Phys. A: Math. Gen. 21, 1173 (1988)), 
an analysis of the smoothing dependence of the correlation functions provides a sensitive means of 
establishing consistency with random matrix theory. We use a semiclassical approach to describe 
these fluctuations and offer a detailed comparison between numerical and analytical calculations for 
an exhaustive set of two-point correlation functions. We consider parametric correlation functions 
with an external Aharonov-Bohm flux as a parameter and discuss two cases, namely broken time- 
reversal invariance and partial breaking of time-reversal invariance. Three types of correlation 
functions are considered: density-of-states, velocity and matrix element correlation functions. For 
the values of smoothing parameter close to the mean level spacing the semiclassical expressions and 
the numerical results agree quite well in the whole range of the magnetic flux. 
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I. INTRODUCTION 

Disordered quantum systems in the metallic regime ex- 
hibit irregular fluctuations of eigenvalues jfl , eigenfunc- 
tions H and also of matrix elements HQ. Parametric 
fluctuations have been discussed in j^]. In the metallic 
regime, which is characterized by a large conductance 
j > 1, such fluctuations can be described by random 
matrix theory (RMT) || on energy scales smaller than 
the Thouless energy E-q = gA, (A is the mean level spac- 
ing). Alternatively, semiclassical methods may be used 
in this regime, as suggested in (j. A semiclassical es- 
timate for parametric correlations of level velocities is 
given in || . Matrix element correlations are discussed in 
P,p|JTo|]. Within a semiclassical approach it is essential 
to incorporate level broadening and work with smoothed 
correlation functions. The level broadening e needs to be 
larger or of the order of the mean level spacing, e > A. 
This ensures that the periodic orbit sums are truncated in 
such a way that only orbits with periods T p shorter than 
the Heisenberg time = 2nh/ A contribute. The results 
reported in j^,|^,|llj predict characteristic dependences on 
the smoothing. As pointed out in ]7j], the smoothing de- 
pendence of the fluctuations provides a sensitive means 
of establishing consistency with RMT. The semiclassical 
approach provides a natural approach of incorporating 
such a smoothing. 

In this paper, we report on extensive numerical cal- 
culations of correlation functions in the two-dimensional 
(2D) Anderson model of localization (l2) in the metallic 
regime, compare also |l3| . In the limit of large g, the sta- 
tistical properties of eigenvalues and eigenvectors in this 
model can be described by RMT on energy scales smaller 
than the Thouless energy Ejy = gA, compare |lj Q. 

We calculate parametric correlation functions (where 
an Aharonov-Bohm flux is used as an external parame- 
ter) as well as fluctuations in systems with weakly broken 
time-reversal (T)-symmetry. We calculate three types 



of correlation functions, namely correlations of the den- 
sity of states Ea|, of velocities [^|JT^1 and of matrix el- 
ements T-invariance is broken by means of 
an Aharonov-Bohm flux (ft. According to RMT, fluctu- 
ations in a T-invariant system, where </> = 0, follow the 
statistics of the Gaussian orthogonal ensemble (GOE). 
At (ft — 0o/4, where </>o = hc/e denotes the flux quan- 
tum, T-invariance is fully broken, and RMT predicts 
the behavior of the Gaussian unitary ensemble (GUE). 
For (ft <C (fto/4: T-invariance is only weakly broken. In 
this case the correlation functions are described by the 
Pandey-Mehta ensemble |l7j . The effect of a weak mag- 
netic field can be exhibited particularly transparently 
within a semiclassical approach. 

All correlation functions calculated in the following will 
be expressed in terms of smoothed spectral densities. In 
the literature, Lorentzian || as well as Gaussian broad- 
ening J?]!! have been used. For numerical calculations, 
Gaussian broadened densities are much more convenient, 
since one invariably deals with finite stretches of spectra, 
and boundary effects are less pronounced due to faster 
decaying tails in the Gaussian case. 

We calculate these correlation functions numerically, 
analyze the smoothing dependence in detail, and deter- 
mine the three non-universal constants, namely the mean 
level spacing A, the conductance g and, in the case of ma- 
trix element correlations, the variance <r 2 ff of off-diagonal 
matrix elements. We report on successes of and problems 
with the semiclassical approach in describing correlations 
in the Anderson model in the metallic regime. 

The article is organized as follows. In Sec. |J, we re- 
call those features of the semiclassical approach that will 
be used in the derivation of the correlation functions. In 



Sec. Ill we describe the Anderson model of localization 
in the weakly disordered regime at finite external flux. In 



Sec. IV we study the correlation functions for the tran- 
sition from the GOE to GUE transition, and in Sec. [V] 
the parametric correlation functions, and compare the 
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study the distribution functions of our results and com- 
pare them to theoretical predictions. We conclude in Sec. 



VII with a discussion of our results. 



semiclassical formulae with the results from the numer- Here, the sum is over periodic orbits p and their repeti- 
ical simulations of the Anderson model. In Sec. VI we tions r. The w PtT are the semiclassical weights, including 

Maslov indices. In general they are complex quantities. 
T p denote the periods and S P (E) the actions of the pe- 
riodic orbits p. Their windings around the flux (f> are 
counted by the winding numbers n p . Similar expressions 
can be derived for densities weighted with level velocities 
or matrix elements as shown e.g. in ^rajO^^j. 

Correlation functions of the type (o) thus involve dou- 
ble sums over periodic orbits. It is argued pl| that the 
average {• • -)e suppresses the non-diagonal contributions 
to this double sum. This is certainly the case for e > A. 
Within the diagonal approximation which amounts to ne- 
glecting the non-diagonal contributions we obtain 



II. THE SEMICLASSICAL APPROACH TO 
UNIVERSAL CORRELATIONS 

In this paper, we calculate correlation functions of the 
following densities. We consider the density of states 
defined as 



d(E, c^) =Y^5 e [E - E a {4>)} 



(1) 



Here, E a (<j)) are the quantum eigenvalues and S e (E) = 
(\/27re)~ 1/2 exp (—E 2 /2e 2 ). Second, we consider the 
density of parametric velocities |lS,n8ln9| 



d v {E, i 



dE a 
dc/> 



Se[E-E c 



(2) 



After unfolding (see the next section), the average level 
velocity is zero. Third, we compute correlation functions 
involving a density of expectation values 



d m {E,4>) = Y,A aa 8 £ [E-E a 



(3) 



with A aa = (ip a (<fi)\A\ip a (<j>)) , where ipaW are the eigen- 
functions corresponding to E a ((f>). A is an operator 
of some real-space observable, not commuting with the 
Hamilton operator H. It is assumed that (A aa ) = 0. 

In all cases, the corresponding densities are decom- 
posed into a smooth and an oscillatory part, 



d{E,<j>) = {d(E,<t>))+d{E,c]>), 



(4) 



where the first term denotes a mean contribution, and 
the second term is a fluctuating part which vanishes upon 
disorder averaging. The mean parts of the densities (j^) 
and (^) are approximately zero. 

For all three densities, we calculate correlation func- 
tions of the type 



E ■ 



(5) 



The average (■ ■ -}e denotes an appropriate average, e.g., 
over disorder realizations and/or energy in the metallic 
regime. Semiclassically, such correlation functions can 
be calculated using a representation of the densities in 
terms of the classical periodic orbits M , 



d(E, 
1 



(G) 



T p exp 



-rSp(E) + 2-nirnp 



s 2 r 2 T 2 



2h 2 



C(4>u<t>2)^j^Y.MT 2 e-^n/^ 



pr 



x I e p *o 
We can then make use of the sum rule [^3| 

J2\^ P \ 2 T^f(Tp)c£dTTf(T). 



(7) 



(8) 



which is valid when long periods T„ dominate the sum 
in (|?]). In order to apply (||) to (0), two further ap- 
proximations are necessary. First, repetitions are ne- 
glected, the usual argument being that periodic or- 
bits proliferate exponentially. Second, assuming that 
the winding numbers are Gaussian distributed, Eq. (0) 
is averaged over the distribution of winding numbers 
P(n,T) = (2ttAT)- 1 /2 eX p(-n 2 /2AT) @. The param- 
eter A = 2T>/L 2 , where V is the diffusion constant and L 
the system size. Evaluating the discrete average over the 
winding numbers by Poisson summation, we then obtain 
the desired semiclassical expressions. 

We remark that the level broadening used in Eq. (0) 
ensures that the periodic orbit sums are truncated in such 
a way that only orbits with periods T p shorter than the 
Heisenberg time in = 2irh/A contribute. We note that 
one could alternatively use a Lorentzian broadening [0. 
For numerical calculations, Gaussian broadened densities 
are much more convenient, since one invariably deals with 
finite stretches of spectra, and boundary effects are less 
pronounced due to faster decaying tails in the Gaussian 



case. 



III. THE 2D ANDERSON MODEL OF 
LOCALIZATION 

We performed numerical simulations within the 2D 



Anderson model of localization 12 1, by diagonalizing the 
Hamiltonian with the help of the Lanczos algorithm [^5| . 
In the site-basis the model Hamiltonian with periodic 
boundary conditions is 



2 



H 



E i n > e »( r 



\n)t r 



(9) 



where \n) represent the Wannier states at sites n in the 
N x N lattice. The on-site potential energies e n are 
taken to be uniformly distributed between — W/2 and 
+W/2. The hopping parameters t nm are non-zero only 
for nearest-neighbor sites n, m and we set the energy scale 
by choosing t = 1 for these sites. For convenience, we as- 
sume that the 2D model is embedded in 3D and defines 
the xy-plane. 

In the presence of a magnetic field the hopping pa- 
rameters acquire an additional factor exp i27r</>/(0o-^O, 
where </> is the magnetic flux, which a periodic orbit en- 
circles in the hopping direction. This phase represents 
the Aharonov-Bohm effect on the system with periodic 
boundary conditions under the magnetic flux. We use 
two magnetic fluxes, <px and 4> y , corresponding to x- 
and y-directions. The corresponding phase of the hop- 
ping parameter is exp i2ir[4> x / '(4>oN) + <fi y /(4>oN)}. For 
completeness, we also study the influence of a homoge- 
neous magnetic field B in z-direction. In this case the 
hopping parameters arc multiplied by exp^fi2itBr y /(j)Q, 
when, e.g., hopping in x-direction. r y is the y-coordinatc 
of the site, and the sign is different for opposite hopping 
directions. To maintain the appropriate periodicity of 
the boundary conditions, B/(j)Q must then be chosen as 
an integer multiple of 1/A^. The hopping parameters in 
y-directions do not change due to B, when we choose the 
vector potential A in the Landau gauge A = (0, Bx, 0). 

The energy spectrum for a single realization of disor- 
der still has an energy dependent density of states. In 
order to study the universal fluctuations, we thus need 
to "unfold" the spectrum |L4| , such that the original set 
of eigenvalues {E a } is mapped to a new set {e Q }, where 

e a = (Af{E a , <j>)) = M(E a , 4>) - ti(E a , 0), (10) 

where Af(E, <ft) — J dE'd(E\ cj>) is the integrated den- 
sity of states, and Af(E a ,(j>) is the fluctuating part of 
Af(E, (j)). In practise we computed (Af(E a , <f))) by fitting 
the Af(E a , 4>) data to a second order polynomial. Then 
we set the value of the polynomial at E a to (J\f(E a , </>)). 
This procedure works particularly well for a relatively 
small number of eigenvalues, where the mean level spac- 
ing A is almost a constant. After unfolding, we have 
A = 1. We remark that an unfolding based on a cu- 
bic spline interpolation |l4| does not work so well in the 
present case. 

The semiclassical approach applies to weakly disor- 
dered systems and for parts of the spectrum, where the 
electron states spread throughout the system. Thus the 
conductance g — tu/tD, with tjj = L 2 /ttT> the Thou- 
less time, should obey j > 1. However, in the infinitely 
large 2D Anderson model, it is well-known that all states 
are localized for any finite amount of disorder pq,E7|. 



Nevertheless, for suitably weak disorder and at small sys- 
tems, one can find large regions in the spectrum for which 
g ^> 1 [ pH , such that we need not go to higher dimen- 
sions to test the semiclassical results. With zero flux 
the (unfolded) spectral fluctuations of the 2D Anderson 
model in this limit of weak disorder are described by the 
GOE of RMT @[l||U. Upon increasing the flux there 
is a transition to GUE ]28| ]. In order to test that we in- 
deed are investigating a part of the spectrum in which 
universality holds, we calculate the nearest-neighbor en- 
ergy level spacing distribution and check that the statis- 
tics for zero flux is given by the Wigner-Dyson result for 
GOE, whereas for finite flux or magnetic field we have 
the GUE result |J . In the following sections we consider 
the dependence of the spectral statistics on the magnetic 



flux 



in the x-direction. The magnetic flux (j> y in 



the y-direction and the homogenous magnetic field in z- 
direction are used as convenient switches between GOE 
(<f> y = and B = 0) and GUE {<j> y ^ or B ^ 0) behav- 
ior. We note that for weak magnetic flux (</>, <j) y <C 0o/4), 
time-reversal symmetry is only weakly broken and the 
statistical properties of the spectrum are described by 
the Pandy-Mehta ensemble |l7p|. 



IV. THE GOE TO GUE TRANSITION 

In this section, we will study the correlation functions 
of the density of states Cd{<j>), the density of level veloci- 
ties C v ((j>), and the density of matrix element correlations 
C m (cj)) as functions of the external magnetic flux = <p x . 
Hence, we also have 4> y = 0, B = 0. We shall always first 
consider the semiclassical derivation of these correlations 
and then turn our attention to a numerical computation 
within the 2D Anderson model. 



A. Density of states 

We first consider correlations of the density of states, 
as defined in Eq. (ry), and calculate the statistic 



Cdfa) = ( d(E, 



(11) 



where (• • -)e denotes an average over a suitably chosen 
energy interval as explained in the last section. Within 
the diagonal approximation H we obtain 



c d (4>) 



- ^zexp(z 2 )erfc(z)} (12) 



with z = [y - 2(f>/(j)o) 2 /5 2 , 5 2 = e/n 2 Xh and erfc(z) the 
complementary error function p9[ . This expression de- 
scribes the crossover of the spectral properties from GOE 
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to GUE behavior, as the flux <f> is varied. A correspond- 
ing expression for a transition driven by a magnetic field 
was given in |3(J. Note that Eq. ([lj) is periodic in (f> 
with period (f>o/2. Eq. ( fl2"| ) can be further simplified in 
the limit of small S (with e > 1). We consider two cases, 
namely cf> — and <f> — </>o/4. In the first case, the sys- 
tem exhibits fluctuations described by the GOE, in the 
second case the fluctuations are described by the GUE. 
We then have 



(13) 



where (3 = 1 in the GOE and (3 = 2 in the GUE. It must 
be emphasized that one requires e > 1 for Eq. ( |l2| ) to 
hold. This ensures that only orbits with periods T p < i H 
contribute to (||). For small values of the level broad- 
ening, the diagonal approximation used in deriving ( |l2] ) 
ceases to be valid pl| . On the other hand, in the limit 
of e <C 1, one has [ B2| 

C{4>) ^(J2*1[E- E a {d?)]) * , (14) 

which is independent of 4>. In summary, one obtains for 
GOE and GUE 



C{4>) 



— t — for e < e c 
2V7re 

li?? for£>£c 



(15) 



Thus the crossover between these two limiting behaviors 
occurs at e c ~ 7r~ 3 / 2 //3. 



Numerical results for the density of states 

We obtained numerical data from the 2D Ander- 
son model for 90 samples of different realizations of 
disorder with W — 2.4, using flux values <fi/<f)o = 
0, 0.007, 0.014, . . . , 0.497. Larger values are not needed 
because of the periodicity of Cd in (j>o/2. There were 
27 x 27 sites in the system. For each disorder real- 
ization we computed 100 subsequent energy eigenvalues 
Ei G [—3.4,-1.9], thereby avoiding contributions from 
localized states in the band tails and from nearly ballis- 
tic states at the band center. We remark that the mean 
density of levels is already nearly constant for this in- 
terval and thus the second order polynomial is ideal for 
the unfolding procedure. After unfolding these eigenval- 
ues, we calculated the Wigner-Dyson statistics P{s) for 
nearest- neighbor level spacings. As shown in Fig. ^, we 
find for flux <f> — that P(s) follows the GOE behavior. 
For flux values close to </>o/4, we have P(s) of the GUE. 
Thus with this choice of parameters we are indeed in the 
ergodic regime of the model as required. 

The comparison between the results for the Ander- 
son model, averaged over all disorder realizations, and 



the semiclassical approximation with different broaden- 
ing values e in units of A is shown in Fig. |^. The agree- 
ment is the best for e > 1, as expected. For smaller 
values there are deviations near the GOE cases <f> = 
and 4> — 0.5</>o- The constant A = 1.21 used in plotting 
Fig. was determined from the statistics of level veloci- 
ties, as we explain below in section IV B. We emphasize 
that in Fig. || and throughout the rest of this paper, we 
have not symmetrized our data with respect to the peri- 
odicity in 4>q . Thus the slight deviations from periodicity 
at 4>q/2 reflect the accuracy of our data. 

In Fig. H we show the small e-behavior of Cd- The 



crossover, predicted in Eq. (|15|) at e c 



-3/2 



0.18//?, can be seen to occur between the values 0.03 < 
e < 0.8 for the GOE, and 0.03 < e < 0.15 for the GUE. 
The upper limits of the intervals in each case can be con- 
sidered as lower boundaries for the validity range of the 
diagonal approximation. The upper validity range of the 
diagonal approximation can also be inferred from Fig. || 
to be close to 4.5 for GOE and 1.7 for GUE. 



B. Density of level velocities 

Next we consider fluctuations of the density of level 
velocities, and compute the statistic 



C v (0) = [ dv 
Within the diagonal approximation, we obtain 



(16) 



V — — oo 
2 



[l + 4z 2 ]V¥exp(z 2 )erfc(z)-4z| (17) 



with z and S as in Eq. (|12|). For small 5 (and with e > 1), 
one obtains the following limiting behaviour 



C v (0) 



for <P = , 
^/irXh/e for = 0q/4 . 



(18) 



Alternatively, in the limit of very small e, we obtain in 
analogy with Eq. ( fl4"|) 



C v (0) 



. (19) 



where fi% ag is the variance of the level velocities dE a /d(j>. 
For (3 = 2, we have @ 



ML K (£)=Moff(£). 



where 




E a ~E a ,~E 



(20) 



(21) 
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With /j,^ s (E) = 2nhX (see section V) we obtain for the 
GUE case (/? = 2) 



C v (0) = 



(22) 



This implies that the semiclassical result of Eq. (P-8|), ob- 
tained within the diagonal approximation, remains valid 
for small e [as opposed to the estimate (fl4|)1. We re- 
mark that while this is true for the GOE (<fi = 0, <fio/2) 
and GUE (0 = 0q/4) cases, it is no longer true in the 
transition regime ]2q] . It will be seen in the next section 
that similar arguments apply to fluctuations of matrix 
elements. 



Numerical results for the density of level velocities 

Using the same data as for the density of states cor- 
relations, we computed C v for the Anderson model with 
different broadenings, as shown in Fig. ^. In this case the 
agreement with the semiclassical approximation is good 
even around = and — 0o/2, i.e. in the GOE case. 
We remark that the shoulders visible in Fig. ^ around 
= 0.0700 and O.430o for the semiclassical expressions 
at small e are an artefact of our approximation for e < A. 

The parameter A was determined from the small e- 
behavior of the C v by fitting the numerical results to Eq. 
(^) as shown in Fig. ||. The agreement of the small 
e-behavior of the numerical data with Eq. ( 5^ ) is rather 
good. Indeed, the agreement is good for all values of e, as 
expected from Eq. (|18j) and discussed above. An alterna- 
tive way is to compute a histogram for the level velocities 
in the unitary case and to use Eq. ( |20| ) and the estimate 
fi^ s (E) — 2ttH\. This is shown in Fig. 0. Both methods 
do not give exactly the same value of A due to numerical 
accuracy and the limited number of samples. With the 
former method we estimate a value A = 1.2 ± 0.1, and 
with the latter one A = 1.4 ± 0.2, where the error limits 
represent the standard deviation of the values obtained 
for different realizations of disorder. We have chosen the 
value A = 1.21 such that the overall agreement of each 
correlation function in Fig. [i] is as good as possible for all 
and all e > 0.3. We emphasize that such an agreement 
is very sensitive on the actual value of A chosen. Fur- 
thermore, we need to assume that A remains constant for 
all 0. As we will show later, this assumption is at least 
questionable for the Anderson model. 



C. Density of matrix elements 

In this section we turn to fluctuations of expectation 
values and consider the statistic 



C m (0) = ( d m (E 



(23) 



and obtain, again in the diagonal approximation, 



c m (<f>) 



E 



2V^e ^ 

v v— — oo 

jexpQ) erfc (^) +exp(z 2 )erfc(z)| (24) 



with z and S as in Eq. (|12]). Moreover, a^ s (E) is the 
variance of non-diagonal matrix elements 



al s {E) = <L4 Qa ,| 2 ) 



(25) 



Correspondingly, <J^ iag (E, 0) is the variance of diagonal 
matrix elements. Unlike it depends on the value of 
the flux 0. In the limiting cases of GOE and GUE, the 
variances are related as 



a 2 di ^(E,<l>) = -a 2 oS (E). (26) 
In the limit of small 6, one obtains for GOE and GUE, 

(27) 



(3 2^e ' 

We shall now argue that these results, derived assuming 
e > 1, remain valid in the limit of small e. Proceeding as 
in the previous section, we obtain for small e 



diag 



2 °$B 



(28) 



2v^Fe (3 2y^e ' 
which is the same as Eq. ( p7j ) calculated for e > 1. 

Numerical results for the density of matrix elements 

We computed eigenvalues and the expectation values 
of the diagonal matrix elements x nn for the dipole mo- 
ment operator x in the site-basis for 69 different realiza- 
tions of disorder W — 2.4 in the Anderson model at flux 
0/^0 = 0, 0.007, . . . , 0.497. We obtained C m with differ- 
ent broadenings e as shown in Fig. Here the agreement 
is reasonable, but not as good as in the two previous 
cases. We note that the small cf> behavior is much bet- 
ter described by the universal v — term than by the 
complete expression of Eq. (|24|) . 

We emphasize that for the present correlation, we 
had to determine two constants describing the system, 
namely, A and a^ s . This makes it even more important to 
have various independent ways of computing them. The 
variance of the off-diagonal matrix elements can be 
determined from the diagonal elements in a similar way 
as the determination of the diffusion constant from the 
level velocities. Namely, we can use the small e-behavior 
of Eq. (27) as shown in Fig. H Interestingly, we find that 
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although Eq. ( p7j ) is expected to remain valid for e ~ 1, 
there are already strong deviations of our numerical data 
from the behavior predicted by Eq. (p7|). This may in- 
dicate that the approximations used in the derivation of 
Eq. (|24|) are less reliable for the matrix element correla- 
tions than for density of states and velocity correlations. 
We can also use the histogram of the diagonal matrix 
elements as shown in Fig. |^. Both methods give slightly 
different values for cr 2 ff in GOE and GUE, whereas we 
assumed in the derivation of the semiclassical formulae 
that a 2 s is independent on the magnetic flux. For GOE 
we obtain a value around <7 2 ff — 0.7 ± 0.1 and for GUE 



a 2 s = 0.8 ± 0.2. Both estimates are compatible within 
the error limits, though. In Fig. ^ we choose a 2 s = 0.65 
in order to get the best overall agreement between Eq. 
( p4| ) and our numerical results. Also, we have again used 
A = 1.21 as an estimate of 2T)/L 2 as in the previous 
sections. 



Keeping in mind the sensitivity of the expressions (12), 
(|l7|), and (|24| ) to the actual values of A and cr 2 ff , we can 
conclude this section by noting that our numerical data 
for the 2D Anderson model in the ergodic regime show 
the main features predicted for the correlations and con- 
vincingly exhibit the GOE to GUE transition. 



V. PARAMETRIC STATISTICS 

In this section, we will study the parametric correlation 
functions of the density of states Kd(A(f>), the density of 
level velocities K v (A(fi), and the density of matrix ele- 
ments K m (A(f>) as functions of the difference in external 
magnetic flux A<p = A(f> Xl averaged over different flux 
values <fi. Since, as studied in the previous section, the 
spectral properties change from GOE to GUE as cf> is 
varied, we introduce an additional flux <f> y = 0o/4 in the 
transverse direction, so as to have spectral statistics ac- 
cording to the GUE for all values of <f>. Again, we shall 
start by first considering the semiclassical derivation of 
these parametric correlations and afterwards compare to 
numerical data from the 2D Anderson model. 



A. Density of states 

For the parametric case we define [jl3| 

K A {A4>) = (d(E, 4>)d*{E, <f> + A<t>)) E ^ , (29) 

where (• • -)e,4> denotes an average over E and (f>. One 
obtains within the diagonal approximation 

oo ^ 

K d {A4>)= J2 7^{l-\/^exp(z 2 )erfc(z)} (30) 



Numerical results for the density of states 

We computed 69 realizations of disorder for the 2D 
system with 27 x 27 sites and a disorder strength W — 
1.7, using the same part of the spectrum as previously 
and flux values cfr/fo = 0, 0.01, 0.02, . . . , 1.0. P(s) reflects 
the GUE, as in Fig. [|, for all values of 4> due to the 
additional transverse flux <p y . In Fig. [To], we show the 
comparison between the semiclassical expression ( |30| ) and 
the numerical data. The agreement is very good for all 
values of e. 

The parameter A was determined in the same way as 
in the GOE to GUE transition in section IV. Because 



the system had been made unitary by introducing an 
additional flux <p y , Eq. ( |22] ) is valid for all values of 4>. 
Consequently, the fitting procedure for the small e- values 
should give the same A for all the flux values, and the 
histogram of the level velocities should have the same 
variance. However, we found differences, which cannot 
be explained only by the error bars. This has been illus- 
trated in Fig. The value A = 2.5, used in Fig. [l^ was 
chosen such that the agreement is the best for all Acf>, all 
e and all three parametric correlations. 

We also used W = 2.4 as in section IV for the GOE 



v— — oo 



to GUE transition and computed the parametric correla- 
tions. But in this case the agreement between the semi- 
classical theory and the data obtained from the Anderson 
model is slightly less convincing than with W = 1.7. 



B. Density of level velocities 

For the parametric correlation of the density of level 
velocities, we define Jl3| 

K V (A4>) = {d v {E,(j>)dl{E,<f> + A<f>)) E ^. (31) 

Within a semiclassical approach, we obtain 

K v (A0) = ^ ± 

V — — 00 

x{(l + 4z 2 ) % A : exp(z 2 )erfc(z) - 4z} (32) 

with z and S as in Eq. (|3^) and for Gaussian broadening. 
This expression is periodic in Ac/) with period </>o- It has 
previously been derived in [|| , using Lorentzian broaden- 
ing, see also . Comparing the v — term of Eq. ( |3^ ) 
with the corresponding expression 

^ffM((i + 4 z 2 )^exp(z 2 )erfc(z)-4z} (33) 

obtained from a Brownian motion model 0], we have 
fil s (E) = 2ttHX (compare section IV). 



with z = {v + A(j)/(j)oy/S 2 and 5 2 = e/ir 2 h\. 
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Numerical results for the density of level velocities 



VI. DISTRIBUTIONS 



Using the same data as in section VA for the density 



of states, we computed the parametric statistics for the 
density of level velocities for the Anderson model. The 
comparison with Eq. ( |32] ) can be seen in Fig. [l^. The 
agreement with the semiclassical approximation is again 
very good. We remark that the overestimation of the 
minima in K v around A<j) = O.10o and O.90o for the 
semiclassical expressions at small e < 0.1 is an artefact 
of the diagonal approximation [Q. 



C. Density of matrix elements 

Lastly, we consider the parametric correlation K m (A<p) 
of matrix elements, i.e., 

K m (A<t>) = (d m {E, (f>)d* m (E, cp + A4>)) E ^ . (34) 

As before we obtain l| 



2 00 

K m (A<t>) = £s. V cxp(z 2 ) erfc(z) 



(35) 



within the diagonal approximation and with z and 5 as 
in Eq. (|30|). We have assumed that the mean density 
of states (d) is essentially energy- and flux- independent. 
Moreover, we have neglected the energy-dependence of 
the off-diagonal variance. 



The distributions of level velocities |jl3[ , shown in 
Fig. |(| and of the diagonal matrix elements of the dipole 
moment operator, in Fig. |9j are well approximated by 
Gaussian distributions, as predicted in random matrix 
theory. According to Eq. the variance of the ma- 
trix elements in the GOE case (<f) — 0) should be ap- 
proximately two times larger than in the GUE case 
(0 w o /4). We obtain a factor of c^agO = 0)/ ( r| iag (^ w 
00/4) « (1.3 ± 0.2)/(0.8 ± 0.1) = 1.6 ± 0.5 in agreement 
with this prediction, although the standard deviations 
are quite large. We again emphasize that the level spac- 
ing distributions obey the Wigner-Dyson statistics, pre- 
dicted in random matrix theory, as shown in Fig. |l| for 
all the disorders and magnetic fields chosen in our work. 

We also calculated the distributions of the off-diagonal 
elements A aa i with E a ~ E a i, We find that their dis- 
tribution is also well approximated by a Gaussian as 
shown in Fig. [l4|. The corresponding variance a^ s should 
be independent of the magnetic flux. This is approxi- 
mately true for our data. With disorder W = 2.4 we get 



'off 



0.8 ± 0.2 in GOE 



0), and 0.9 ±0.2 



in GUE 
find 



'off 



0o/4, (f)y = 0) and with W = 1.7, we 
= o /4 and 0.6 ± 0.2 



= 0.7 ± 0.2 at 6 T = 0, 



at cf> x = <j)y = 0o/4. The error bars represent again the 
standard deviations of the values obtained for different 
realizations of disorder. 



Numerical results for the density of matrix elements 

In Fig. [l3L we show the comparison between semiclas- 
sical and numerical results for the parametric statistics of 
the matrix elements of the dipole moment operator using 
the same data as for the two previous parametric correla- 
tions. The agreement here is even better than in the GOE 
to GUE transition. This is noteworthy, because of the 
large discrepancies between the values of A for different 
flux values (cp. Fig. ^) which we neglected in the semi- 
classical derivation of Eq. ([55]). The off-diagonal variance 
a^g was determined in the same way as in section ^ for 
the GOE to GUE transition, giving a% s = 0.50 ± 0.05. 
We get different values for different flux values as for the 
diffusion constant, but the variations are much smaller. 
By calculating directly the variance of the matrix ele- 
ments between nearest-neighbor sites we get a slightly 
larger value c 2 ff — 0.65 ± 0.05. Here, the error bars rep- 
resent the deviations from the average value for different 
flux values. 

Thus in summary, we find that as in section the 
general behavior of the data obtained for the 2D Ander- 
son model in the GUE case is very well reproduced by 
the semiclassical expressions (|30|), (B2J) , and (35). In fact, 



the agreement is even better than in section [V 



VII. CONCLUSIONS 

In this paper, we have reported on extensive calcula- 
tions of smoothed correlation functions in the 2D Ander- 
son model of localization. We have calculated correlation 
functions of energy levels, their parametric derivatives 
and of diagonal matrix elements in the metallic regime 
(g 3> 1). For two cases, namely for parametric corre- 
lations and for fluctuations in the transition regime be- 
tween GOE and GUE, we have presented detailed com- 
parisons of our numerical results with semiclassical the- 
ory, focussing on the dependence of the fluctuations on 
the level broadening. 

Our results can be summarized as follows. First, one 
expects the semiclassical theory to be appropriate for 
level-broadenings in the range of 1 < e -C g (with e in 
units of A). Comparison with asymptotic expressions for 
small e [Eqs. (§§) and @] shows that the lower 

bound actually extends to e c ~ 7r~ 3 / 2 //3 for density-of- 
states fluctuations. In the case of fluctuations of level 
velocities and matrix elements, moreover, the diagonal 
approximation remains valid for arbitrarily small e. This 
is simply due to the fact that the additional factors in 
Eqs. (Q) and (^) are essentially random and help to sup- 
press off-diagonal contributions to (|?j). Our numerical 
results verify these conclusions. 
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Second, at large values of e we observe deviations from 
the universal theoretical results, as expected. This is evi- 
dent in Figs. ||, [| and |[ The value of the conductance in 
this case is g — 12 ± 3. Interestingly, in Fig. || in partic- 
ular, we observe deviations from the universal prediction 
at considerably smaller values of e. From this we conclude 
that fluctuations of matrix elements are particularly sen- 
sitive to non-universal effects. This is consistent with 
the following observation. In the universal regime, the 
semiclassical expressions derived in this paper should be 
dominated by those terms for which \ v — 20/0o| is mini- 
mal. However, in the case of matrix element fluctuations, 
non-universal contributions are particularly large (com- 
pare Fig. Q). This is not surprising since it can be shown 
that short periodic orbits make large, non-universal con- 
tributions to C m (</>). 

Third, in the case of parametric fluctuations (Figs. [l(], 
|l2| and [l|) we observe excellent agreement with the semi- 
classical predictions. This is due to the fact that (i) 
these numerical results are averaged over a considerably 
larger ensemble and (ii) that the conductance is larger 
(g = 24 ±5). 

Fourth, we emphasize that in our case the parameters 
g and <t^ s (E) are found to depend on the magnetic flux 
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derforschungsbereich 393. 



(compare Fig. 11). The flux dependence turned out to 



be more prominent with the smaller disorder strength 
we used. That is why our numerical results for the cor- 
relations in the GOE to GUE transition agrees better 
with the semiclassical formulae with W — 2.4 than with 
W = 1.7, even if the conductance is smaller in the former 
case. Within the framework of the semiclassical theory g 
and <Tq S (E) are expected to be independent of <fi since an 
Aharonov-Bohm flux does not change the classical me- 
chanics. 

Fifth, we have verified the relation between the vari- 
ances of diagonal and non-diagonal matrix elements in 
the GOE and GUE. The agreement of our numerical re- 
sults with the prediction is reasonably good . 

In summary, we have shown to which extent fluctua- 
tions in the 2D Anderson model are accurately described 
by universal semiclassical formulae. We have found, in 
particular, that the fluctuations depend sensitively on 
the level-broadening and that this dependence can be 
used to assess consistency with RMT, as originally sug- 
gested in M. This is particularly important for the fol- 
lowing reason. In order to test recent predictions [j3Clf on 
the effect of incipient localization on the fluctuations of 
wave-function amplitudes in the 2D Anderson model it 
is essential to have an accurate and quantitative under- 
standing of the metallic regime. 
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12 3 



FIG. 1. Histograms for energy level spacings of the un- 
folded energies for all samples with disorder W = 2.4 and a 
system size TV 2 = 27 2 . The (smooth) lines denote the GOE 
(solid) and GUE (dotted) Wigner-Dyson distributions || for 
= and (f>/<f)o = 0.25 ± 0.05, respectively. 
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FIG. 2. GOE to GUE transition for density of states cor- 
relations according to Eq. ( |l2j ) (solid lines) and corresponding 
results from the numerical simulations of the Anderson model 
(symbols). The parameter A = 1.21 and e = 0.316 (<), 0.447 
(o), 0.631 (□), 0.891 (O), 1.26 (A), 1.78 (+), 2.51 (*). 
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FIG. 3. Small e-behavior of C d with W = 2.4. The solid 
line indicates l/(2-7r 2 e 2 ), the long dashed line is l/(47r 2 e 2 ) and 
the short dashed line denotes l/(2y / 7re). 
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FIG. 4. GOE to GUE transition for level velocity corre- 
lations according to Eq. ( |l7j ) (solid lines) and corresponding 
results from the numerical simulations of the Anderson model 
(symbols) . The parameters are the same as in Fig. ^. We ad- 
ditionally include the broadenings e — 0.158 (v) and 0.224 
(x). 



1CT 2 1CT 1 10° 

8/A 

FIG. 5. Determination of the parameter A from the 
small e-behavior of G v at W = 2.4. By fitting the four 
first points (■) on the left to Eq. (^) one gets a value 
A = 1.22±0.02, whereas fitting the first eight points (bold □) 
gives A = 1.30 ± 0.01. The difference between a plot of Eq. 
@ with A = 1.22 (thin dashed line) and A = 1.3 (thin solid 
line) is very small. 



0.15 




8E,/&|> 

FIG. 6. Distribution of level velocities averaged over flux 
values (f>/<j>o = 0.175, . . . , 0.329 and 90 different realizations of 
disorder for W = 2.4. The line represents a fit by a Gaussian 
distribution. 
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0.2 0.3 0.4 0.5 
<H>0 

FIG. 7. GOE to GUE transition for matrix element corre- 
lations according to Eq. ( p^ ) (solid lines) and corresponding 
results from the numerical simulations of the Anderson model 
(symbols). The parameters are the same as in Fig. ^ and the 
off-diagonal variance is taken to be a^ ff = 0.65. The dashed 
lines indicate the v = term of Eq. (E4) for small <p. 




FIG. 8. Determination of the off-diagonal variance of the 
dipole moment operator by fitting the values of C m at small 
e (filled symbols) to Eq. (pi). The disorder is W = 2.4. 




FIG. 9. Distribution for diagonal matrix elements of 
the dipole moment operator in units of the lattice constant 
and in case W = 2.4, = (filled circles) and averaged 



over all 



i) 



0.175, ...,0.329 (open circles) with 



0. 



The lines are fits by Gaussian distributions. The variance 
a| iag = 2al s /f3 of the data is 1.3 ±0.2 for cj> = and 0.8 ±0.1 
for (j>/4>o = 0.175, 0.329. The error limits represent the 
standard deviations of the values obtained for different real- 
izations of disorder. 




0.4 0.6 0.8 1.0 

FIG. 10. Parametric correlations of density of states ac- 
cording to Eq. (^) (solid lines) compared to the numerical 
results for the Anderson model (symbols) as a function of Acj). 
The parameters are W = 1.7, A = 2.5 and e = 0.112 (>), 
0.158 (v), 0.224 (x), 0.316 (<), 0.447 (o), 0.631 (□), 0.891 
(O), 1.26 (A). The curves have been shifted by multiples of 
0.1 for clarity. 
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0.0 0.1 0.2 0.3 0.4 0.5 

4»/4» 

FIG. 11. Parameter A = 2V/L 2 for the Anderson model, 
determined by fitting Eq. ( p2[ ) to the data (O) and from the 
variance of the level velocities (■) with W = 1.7 and different 
flux values in the presence of a transversal flux <j> y = </>o/4. 
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FIG. 12. Parametric correlations of level velocities accord- 
ing to Eq. ( |32[ ) (solid lines) compared to the numerical results 
for the Anderson model (symbols) as a function of Acj>. The 
parameters are the same as in Fig. [l(]. The curves have been 
shifted by multiples of 1 for clarity. 
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AW 

FIG. 13. Parametric correlations of matrix elements ac- 
cording to Eq. (solid lines) compared to the numerical 
results for the Anderson model (symbols) as a function of Acj>. 
The parameters are the same as in Fig. [j^ and cr^g = 0.48 
has been used. The curves have been shifted by multiples of 
0.1 for clarity. 
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FIG. 14. Distribution of real and imaginary parts of the 
off-diagonal dipole matrix elements at flux values = (filled 
symbols) and (j> = 0o/4 (open symbols). Additionally, </>y — 
for W = 2.4 (o) and <j) y = o /4 for W = 1.7 (o). The lines 
represent fits by Gaussians. The distributions for W — 1.7 
have been shifted by 0.2 for clarity. 



12 



